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We use first-principle density functional theory (DFT) to study the transport properties of single 
and double barrier heterostructures realized by stacking multilayer h-BN or BC2N, and graphene 
films between graphite leads. The heterostructures are lattice matched. The considered single barrier 
systems consist of layers of up to five h-BN or BC2N monoatomic layers (Bernal stacking) between 
graphite electrodes. The transmission probability of an h-BN barrier exhibits two unusual behaviors: 
it is very low also in a classically allowed energy region, due to a crystal momentum mismatch 
between states in graphite and in the dielectric layer, and it is only weakly dependent on energy in 
the h-BN gap, because the imaginary part of the crystal momentum of h-BN is almost independent 
of energy. The double barrier structures consist of h-BN films separated by up to three graphene 
layers. We show that already five layers of h-BN strongly suppress the transmission between graphite 
leads, and that resonant tunneling cannot be observed because the energy dispersion relation cannot 
be decoupled in a vertical and a transversal component. 

PACS numbers: 81.05.ue,72.80.Vp,73.63.-b,73.50.-h,71.15.Mb 



Mobility of suspended graphene can be extremely high, 
as demonstrated by experiments and theor}*^ — . However, 
use of graphene in solid state devices typically requires 
deposition or growth on a dielectric substrate, which can 
strongly suppress mobility, due to electron coupling with 
dielectric-layer phonon modes, as demonstrated in the 
case of Si02 or Hf02, the most commonly used dielectric 
materials^'^. 

Hexagonal boron nitride (h-BN) has been recently in- 
vestigated as a promising dielectric for graphene^^^, 
while boron nitride domains has been suggested as a way 
to engineer graphene nanoribbon transport propertied. 

h-BN has hexagonal geometry and a lattice constant 
closely matching that of graphene^^, and the electronic 
interaction of h-BN with a single graphite sheetf^i^ or 
with bilayer graphene^^ is very weak. As a result, 
graphene deposited on an h-BN substrate maintains its 
electronic and transport properties. 

From this perspective, understanding transport along 
stacked graphene/h-BN structure can provide relevant 
information regarding the leakage current in graphene 
based Field Effect Transistors exploiting h-BN as gate 
dielectric or the performance of recently proposed vertical 
devices^^^^^. 

In this work, with ab-initio method, we calculate the 
quantum transport in single and double barriers con- 
sisting of h-BN and BCN layers between semi infinite 
graphite leads, using the formalism introduced by Choi 
et al.^^. Barrier consists of up to five h-BN or h-BC2N 
atomic layers arranged as shown in Fig. [TJ according to 
a recent paper by Ribeirc^^. The double barrier is made 
of two or four h-BN or h-BC2N layers separated by up to 
three graphene layers. The scattering region considered 
in the calculation includes at least four atomic layers of 
graphite on each side of the barrier film, to ensure that 
the charge density does not change at the interface be- 
tween the lead and the scattering region. The geometry 
of the total structure in the scattering zone has been op- 
timized. 
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FIG. 1. Side view of the supercell used to rep- 
resent the scattering region corresponding to a) 
graphite 1 5 (BN) I graphite; b) graphite 1 5 (BC2N)| graphite; 
c) graphite|2(BN)|graphene|2(BN)|graphite. 



Ab-initio calculations have been performed by means 
of Quantum Espresso using a plane wave basis set in 
the local density approximation (LDA)^^. A 35 Ry wave 
function cutoff has been considered, the Brillouin zone 
has been sampled using a 30 x 30 x 30 Monkhorst— Pack 
grid. Atomic positions have been relaxed using a conju- 
gate gradient algorithm until all components of all forces 
was smaller than 0.01 eV/Awhile the electronic mini- 
mization has been performed with a tolerance of 10~^ 
eV. It has been shown^^ that LDA provides a reliable 
description of the geometry and the electronic struc- 
ture in the presence of weak interactions between h-BN 
layers or between h-BN and graphene layers. On the 
other hand, the LDA representation for the exchange- 
correlation potential cannot correctly describe the ex- 
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cited states in organic systems and leads to overestima- 
tion of conductance. However, following other works on 
this subject^^, we expect that these effects have only mi- 
nor consequences, and thus our conclusions should not 
be affected by this choice. The transmission probabili- 
ties have been calculated with the PWCOND^^ module 
of Quantum Espresso. The transport properties are stud- 
ied in the framework of Landauer formalism, ^^^^^ where 
the ballistic conductance is given by G = GqT{E)^ where 
Go is the quantum conductance Go = 2e^//i, and T {E) 
is the total transmission at the energy E. 
T{E) is obtained as 

nE) = Y.J2J2T^Ah,E) (1) 

k\\ i j 
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where Tij (/cy , E) is the probability that an electron with 
energy E and transversal momentum k\\ incoming from 
the i—th Bloch state is transmitted to the outgoing j— th 
state of the other electrode. Sums runs on both spins. 
The first sum is performed over /cy belonging to the two- 
dimensional Brillouin zone (2D-BZ) of the supercell. 

The transmission probability of a single h-BN barrier 
is shown in Fig. [2] for a different number of layers. For a 
given number n of single h-BN layers, the transmission 
probability is almost independent of energy, in an en- 
ergy range of about 3 eV starting from the Fermi level. 
This behavior can be explained by the structure of the 
complex bands in bulk h-BN shown in Fig. [3l in the en- 
ergy interval considered, in the symmetry point K on the 
plane parallel to the layer interface (where the graphite 
states with real momentum reside), h-BN has a disper- 
sion relationship according to which the imaginary wave 
vector is almost independent of energy, and therefore T 
as well. Of course T decays exponentially with increasing 
n, as can be verified by plotting the average value (f ) 
in the energy range between the Fermi energy value and 
3.4 eV above as a function of the number of layers, as in 
the inset of Fig. [2l The fitting exponential is 

f ^ e-^-^^ (2) 



FIG. 2. Top: Density of states of bulk h-BN (solid black line) 
and graphite (dashed red line). Bottom: Tunneling proba- 
bility as a function of the incident energy for a single h-BN 
barrier of n atomic layers between graphite leads 
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FIG. 3. Complex band structure for h-BN bulk. 



and five h-BN layers strongly suppress the transmis- 
sion. The dip in correspondence of £^ — £^Fermi = eV 
is due to the null density of states in graphite leads. In 
addition, we also note that for E — £^Fermi < —3 eV, 
the transmission probability is suppressed even in a clas- 
sically allowed region, probably because of a mismatch 
between propagating states in the h-BN barrier and in 
graphite. 

In Fig. m the transmission probability as a function of 
energy is shown for a single barrier of BC2N of n atomic 
layers between graphite leads. The atomic structure is 
shown in the inset. A qualitatively similar behavior as 
h-BN is observed, with an exponential decay of the trans- 
mission probability in a narrow range around Fermi en- 
ergy. BC2N has a gap of 1.6 eV. 

Fig. [5] shows the transmission probability as a function 
of energy for some single and double barrier structures. 
All these systems are symmetrical, with two identical 



barriers separated by one, two or three graphene layers. 
For clarity, hereinafter, we have indicated the single bar- 
rier systems with the acronym SB followed by the number 
of h-BN layer (e.g. in Fig. [1^, SB5). In the same way, 
the double barrier systems are related with DB and three 
numbers, namely the number of h-BN layers on the left, 
the number of graphene layers in the central region, and 
the number of h-BN layers on the right (e.g. in Fig. [Th, 
DB2|1|2). By comparing the transmission of the single 
h-BN barrier (SBl) with that of the double h-BN barrier 
(DB1|1|1) (Fig. [5^), the noteworthy aspect is that the in- 
sertion of a single graphene sheet has a noticeable effect 
on transport even in the direction perpendicular to the 
plane (Fig. [5j). While the insertion of even a single BN 
(SBl) sheet strongly suppresses the transmittance (par- 
ticularly in the zone near the Fermi energy), the presence 
of a triple layer (DB1|1|1) yields almost the same behav- 
ior as the single-BN sheet . In addition, DB1|1|1 has a 
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FIG. 4. Transmission probability as a function of energy for 
a single BC2N barrier of n atomic layers between graphite 
leads. The inset shows the atomic structure of the BC2N)n 
isomer studied. 



FIG. 6. a) Transmission probability as a function of the en- 
ergy (eV) for DB2|1|2; b) and c) Transmission probability as 
a function of energy for DB2|1|2 system with step sampling 
2 X 10"^ and 4 x 10"^ 
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FIG. 5. Transmission probability as a function of the energy 
(eV) for different systems; a: SBl (solid black line), SB2 
(dashed red line), DB1|1|1 (dotted green line), b: DB1|1|1 
(solid black line), DB1|2|1 (dashed red line), DB1|3|1 (dotted 
green line) 



much higher transmittance than SB2, despite the iden- 
tical number of h-BN layers. In Fig. [Sb, we report the 
transmittance of double-barrier systems in which two h- 
BN layers are separated by one, two or three graphene 
layers. 

It is also important to notice that no resonant tun- 
neling occurs in the double barrier structures, in the 
whole energy range explored, as we have verified with 
a very dense energy sampling of transmission probabil- 
ity, down to steps of 10~^ eV, as shown in Fig. preferred 
to (DB2|1|2) case. The reason is due to the fact that the 
energy dispersion relation in graphene cannot be writ- 
ten as the sum of a longitudinal and a transversal com- 
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FIG. 7. Transmission probability of the double barrier sys- 
tem DB2|1|2 as a function of energy for selected A:||: a) 
(0.28000;0.28867); b) (0.31000;0.32909); c) (0.28000;0.33486); 
in unit — with a=2.503 A 



ponent. Therefore, for each transversal wave vector /cy, 
one observes resonances in T{E) with unity transmission 
probability peaks, as shown in Fig. [71 at different ener- 
gies. When the sum over k\\ resonance do not occur at 
the same energy and are therefore washed out. 



I. CONCLUSION 

In this paper we have investigated the tunneling proba- 
bility of vertical heterostructures consisting of single and 
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double layers of h-BN and BC2N separated by graphene 
sheets and connected to graphite leads. These structures 
are interesting from the point of view of graphene elec- 
tronics, also because few devices have been proposed in 
which transport occurs in the direction perpendicular to 
graphene plane o^^i^^i^^ . As expected, the transmission 
probability is exponentially dependent on the number 
of layers, and is already strongly suppressed by a sin- 
gle monolayer of h-BN or BC2N. In addition, due to 
the energy dispersion relationship in the energy gap, for 
which the imaginary wave vector is almost constant as a 



function of energy, the transmission probability has only 
small dependence on energy. We have also shown that 
resonant tunneling is completely suppressed by the pe- 
cuHar energy dispersion relation of graphene, that can- 
not be decomposed in the sum of a longitudinal and a 
transversal component. 
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